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Abstract. We report the detection of WC stars in five Wolf- 
Rayet (W-R) galaxies: He 2-10, NGC 3049, NGC 3125, NGC 
5253 and Tol 89. The faint broad C IV A5808 Une requires suf- 
ficiently high S/N (> 40) to be detected explaining the non- 
detection of this WC feature in previous observations. From the 
measurement of W-R emission Unes (N III A4640H-C III A4650, 
He II A4686, and C IV A5808), we conclude that all W-R re- 
gions contain a mixed population of WNL, and early WC stars. 
The exception is the high-metallicity region NGC 3049 where 
late WC stars prevail. 

A spatial offset between the multiple peaks of the nebu- 
lar emission and the stellar light in He 2-10 and Tol 89 is 
observed. These nebular emission structures are likely due to 
the existence of bubbles and loops, owing to the injection of 
mechanical energy in the ISM through the W-R winds and/or 
supernovae. Due to age differences and likely smaller energy 
deposition the structures around the W-R regions are possibly 
smaller than the ones predominantly energized by SNe. The 
spatial distribution of W-R stars closely follows the stellar con- 
tinuum with no significant distinction between WN and WC 
stars. 

From the luminosity of the W-R signatures we have esti- 
mated the absolute number of W-R stars of the different sub- 
types. The WCAVN number ratios have typical values between 
0.2 - 0.4, and show no clear trend with metallicity. For low- 
metallicity objects (Z ~ 1/5 Z©), these values are larger than 
the observed WC/WN ratios in Local Group objects, but are 
compatible with expectations for star forming events with short 
duration if stellar evolution models with high mass loss are 
used. 

We derive ages for the starburst regions in the range of 3 
to 6 Myr and confirm that the burst duration must not exceed 
^ 2-4 Myr to account for the high population of W-R stars 
observed in starburst regions, even if emission line stars simi- 
lar to those observed in R136 and NGC 3603 are common in 
starbursts. Within the uncertainties the majority of the observed 



Send offprint requests to: D. Schaerer, schaerer@obs-mip.fr 

* Based on observations obtained at the European Southern Obser- 
vatory (ESO), Chile 

** Present address: European Southern Observatory, Karl- 



quantities is reasonably reproduced by models with a Salpeter 
IMF. Although some W-R lines in few regions are stronger than 
predicted by the models no clear case requiring a significantly 
flatter IMF is found. IMF slopes much steeper than Salpeter 
may, however, not be compatible with our data. 

Key words: galaxies: individual: He 2-10; NGC 3049; NGC 
3 125; NGC 5253; Tol 89 - galaxies: starburst - galaxies: stellar 
content - Hll regions - stars: Wolf-Rayet 



1. Introduction 

Starbursts play a major role in the global process of galaxy for- 
mation and evolution. However, despite two decades of intense 
research on starburst galaxies, our present knowledge of the in- 
trinsic starburst properties is still rather poor and fundamental 
questions remain to be answered. For example, what is the du- 
ration of starbursts ? Which type of stellar population is formed 
during these events ? Is the propoition of massive stars higher 
in starbursts than in more "quiet" star-forming regions? Obser- 
vationaly, answering these questions is not an easy task. In this 
respect the so-called "Wolf-Rayet (W-R) galaxies" may be the 
ideal laboratories since these objects harbour the most massive 
stars known, O stars and their descendents W-R stars, which 
allow to probe the upper part of the initial mass function (IMF) 
and the youngest stellar populations. 

W-R galaxies are usually defined as "those galaxies in 
whose integrated spectra a broad emission feature at Hell 



A4686 attributed to W-R stars has been detected" (Conti, 1991 ) 



In practice, the detection of a broad emission feature at A4650— 
4690 (the so-called "W-R bump", which possibly includes ad- 
ditional emission lines, cf. below) attributed to W-R stars is 
often simply used. 



Whereas the initial compilation of Conti ( 1991 ) includes 37 
galaxies, to date at least 100 such objects are known (see the 
compilation of Schaerer & Contini, 1998). W-R galaxies are 
found among emission-line galaxies. However, the class of W- 
R galaxies encompasses a wide range of properties and galaxy 
types. While most of them fall in the category of Hll galaxies, 
broad W-R emission lines have also been detected in more "ex- 



2 



Daniel Schaerer et al.: Populations of WC and WN stars in Wolf-Rayet galaxies 



1988), Seyfert 2 galaxies (e.g. Heckman et al., 1997) and giant 
cD galaxies, lo cated in the c entre of galaxy clusters with strong 
cooling flows ( 411en, 1995 ). 

The mere presence of W-R stars in these objects indicates 
recent (< 10 Myr) star formation and the existence of massive 
stars (A/initial ^ 25 Mq, cf. Maeder & Conti, 1994); this pro- 
vides already interesting information about star formation in 
these objects. Quantitative analysis of the He II A4686 emis- 
sion or the W-R bump show that star formation must have 
occurred over short timescales compared to the lifetime of 



NGC 5253, Tol 89), three of them ai-e in common with VC92. 
A first account of our observations of NGC 5253 was given in 
Schaerer et al. (|1997|). 



massive stars (|Kunth & Sargent, 198 1^ [Arnault et al., 1989 



Vacca & Conti, 1992| ; [Meynet, 1995[ ), and have allowed to de- 
rive the absolute number of W-R stars present in these regions 
(iVw-R 100 - 10^ e.g. Vacca & Conti 1992, hereafter 
VC92). 

This blue W-R bump is often blended with nearby nebular 
emission lines of He, Fe, or Ar, and can show several broad 
stellar- emission components (N III A4640, C III A4650, He II 
A4686) difficult to separate in most low or medium-resolution 
spectra. Their origin can a priori be due to W-R stars of WN 
and/or WC subtypes, which has been a question of some earlier 
dispute (cf. Osterbrock & Cohen, 1982; Sargent & Filippenko, 
1991; with Kunth & Schild, 1986; Conti, 1991). 

The strongest emission line of WC stars is CiV A5808 
which is very weak in WN stars. This "red W-R bump" 



has quite rarely been observed so far. Although Conti (1991) 



claimed the absence of convincing evidence for carbon features 
in W-R galaxies (cf. also Kunth & Schild, 1986) there is no 
doubt about their existence anymore (e.g. NGC 4861: Diner- 
stein & Shields, 1986; NGC 4214: Sargent & Fillipenko, 1991, 
Kobulniky & Skillman, 1996; NGC 2363: Gonzalez-Delgado 
et al., 1994; Mrk 996: Thuan, Izotov & Lipovetsky, 1996; Mrk 
475, Mrk 1450: Izotov et al., 1994; NGC 7714: Garcia- Vargas 
et al., 1997; I Zw 18: Izotov et al., 1997, Legrand et al., 1997b). 
Where the data is available, C IV A5808 is generally weaker 
than He II A4686. In our recent compilation (Schaerer & Con- 
tini, 1998) we find ^ 15 objects with a fairly well established 
detection of broad C IV A5808. 

The presence of WC stars, although possibly less numerous 
than WN stars, is indeed expected both from observations of 
W-R populations in the Local Group (e.g. Massey & Johnson, 
1998) and from stellar evolution models (Maeder & Meynet, 
1994). The predictions are that, depending on the evolution- 
ary model and metallicity, — 55 % of young starburst should 
show an important WC population (Meynet, 1995; Schaerer & 
Vacca, 1998; hereafter SV98). Most surprisingly the homoge- 
nenuous, fairly high signal-to-noise spectra from the sample of 
VC92 containing 12 regions with He II A4686 detections (and 
10 upper limits) revealed only one C IV A5808 detection (in He 
2-10 A). If true, this would certainly contradict the predictions 
from synthesis models (Meynet, 1995; SV98) using high mass 
loss stellar evolution models, which otherwise compare well to 
observations in the Local Group (Maeder & Meynet, 1994). 

To verify if this apparent discrepancy really holds we have 
initiated a search for WC stars in well-known W-R galaxies. 



In the present paper we report the successful discovery of 
both WN and WC stars in all of the objects, which supports 
our initial hypothesis of an observational bias (i.e. a too low 
signal-to-noise ratio around 5800 A) against their detection in 
the VC92 sample, and (at least partly) relieves the aforemen- 
tioned discrepancy. We then use both the WN and WC sig- 
natures to constrain the burst properties and stellar evolution 
models. 

The paper is structured as follows. The observations and 
reductions are described in Sect. 2. A spatial analysis of the 
emission lines is shown in Sect. 3. We describe the properties 
of the W-R regions and their massive star content in Sect. 4. 
Constraints on the evolutionary tracks and the properties of the 
starburst regions (age, burst duration, IMF) are derived in Sect. 
5 from a comparison with evolutionary synthesis models. Fi- 
nally, our main results are summarised and discussed in Sect. 
6. 

2. Observations and data reduction 

Long-slit spectra of galaxies were obtained on the nights of 
1995 April 24 - 26 at the ESO 2.2m telescope. The data were 
acquired with the EFOSC2 spectrograph and a 1024 x 1024 
Thomson CCD with a pixel size of 0.34". We used grism 4 
which gives a spectral coverage of 4400 - 6500 A with a reso- 
lution of ~ 5 A. During the nights, we also observed the spec- 
trophotometric standard stars HD 84937 and Kopff 27 in order 
to flux calibrate the spectra of the galaxies. Note that only the 
first night was photometric. Spectra of He-Ar calibration lamp 
were obtained immediately before and after the galaxy integra- 
tions in order to accurately calibrate the wavelength scale. 

The slit was oriented in order to cover regions with previ- 
ous detections of W-R stars and others bright optical regions 
in galaxies. The slit width was 1.6" for the galaxy observations 
and 5" for the standard stars. The total integration time for each 
galaxy is given in Table |l]. We took multiple exposures of 1200 
s in duration, short enough to avoid saturation of the bright neb- 
ular lines (H/3 and [O III]) and to recognize cosmic ray impacts. 
The seeing was relatively stable during the observations with 
a mean spatial resolution of about 1", as measured by the La 
Silla seeing monitor. The spectra were acquired at low airmass 
(< 1.1, except for NGC 3049 observed with an airmass 1.3) 
and no correction for the loss of blue light due to atmospheric 
dispersion was made. 

The spectra were reduced according to standard reduction 
procedures using the MIDAS package LONG. A bias level was 
subtracted from each frame using the overscan region of the 
CCD chip. A combination of flat field images from the twilight 
sky and a continuum quartz lamp were divided into each frame 
to coiTect for pixel to pixel variations along spatial and disper- 
sion axes respectively. After wavelength calibration, transfor- 
mation to equal wavelength intervals along the dispersion axis 
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Galaxy 




Coordinates (J2000) 


Type 


rriB 


E(B-V) 


^0 


D 




Exposure time 


P.A. 


Scale 






a 


<5 




[mag] 


[mag] 


[km s~i] 


[Mpc] 


[mag] 


[min] 


n 


[pc/"] 


He 2-10 


OS** 


36" 16=0 


-26°24'40" 


10 Pec 


12.5 


0.21 


873±10 


8.8 


-17.3 


80 (4x20) 


94 


43 


NGC 3049 


09'^ 


54" 49=9 


+09° 16' 19" 


SB(rs)ab 


13.0 


0.04 


1494±4 


18.3 


-18.3 


80 (4x20) 


35 


89 


NGC3125 


lO** 


06" 34M 


-29°56'10" 


E? 


13.3 


0.25 


1080±47 


11.5 


-17.0 


100(5x20) 


123 


56 


NGC 5253 


jjh 


39" 55=8 


-3I°38'4l" 


E/SO? 


11.1 


0.05 


404±4 


4.1 


-17.2 


90 (4x20 -^ 10) 


24 


19 


Tol89 


14h 


01" 25M 


-33°04'28" 


GHlI in SBdm 


16.0 


0.23 


1226±11 


14.7 


-14.8 


120(6x20) 


39 


71 



Global parameters come from RC3 except for the heliocentric radial velocity (Vq) of He 2-10 which comes from Kobulnicky et al. ( 1995 ). The 
Galactic foreground extinction E(B-V) comes from Burstein & Heiles (1984). The distance (D) is derived from the Galactic Standard of Rest 
(GSR) velocity, using a Hubble constant Ho — 75 km s^^ Mpc^^, except for the distance of NGC 5253 taken from Saha et al. (1995). The 
absolute blue magnitude (Mabs) is derived from the apparent magnitude (me) and the adopted distance (D). The position angle (RA.) of the 
slit is measured from North to East. The linear scale (in pc/") is computed using the adopted distance adopted distance (D) 



from the same galaxy were combined to average frames and re- 
ject cosmic ray events. Care was taken to discriminate between 
cosmic rays and the high peaks of strong emission lines which 
might also be confused with cosmic ray events. Tracking with 
the autoguider was sufficiently good that no shifting or registra- 
tion of the frames was required prior to the combine procedure. 
Emission-lines from the night sky were then subtracted using a 
linear fit to at least 50 "sky" pixels on either side of the galaxy. 

3. Detailed spatial analysis 

We first use our long-slit spectroscopic observations to com- 
pare, for each galaxy, the spatial distributions of nebular and 
W-R emission-line intensities and the stellar continuum emis- 
sion. For this purpose, we extracted one-dimensional spectra 
along the slit by adding 3 pixels to match the seeing of our ob- 
servations (^ 1"). We choose one pixel of overlap between two 
adjacent spectra to maximize signal-to-noise ratio while still 
maintaining a good spatial resolution. Depending on the spa- 
tial extent we performed between 80 and 180 extractions per 
galaxy, such that the brightest emission-lines (H/3 and [O III] 
A5007) were strong enough to be measured with reliability. 

We then measured, for each individual spectrum, the pa- 
rameters (central wavelength, FWHM, flux and equivalent 
width) of the brightest nebular emission lines (He I A4471, H/3, 
[Oiii] A4959, [Oiii] A5007, Hex A5876 and [Oi] A6300), the 
broad emission lines due to W-R stars (blue W-R bump: blend 
of N III A4640 and C III A4650, He II A4686; red W-R bump: 
C IV A5808) together with the continuum flux under each emis- 
sion line. 

3.1. Variations of emission line and continuum intensities 

The results of the spatial analysis are displayed in Figs. [l| to 
||. We do not observe any significant variation in the spatial 
distribution of the bright nebular emission lines. We therefore 
plot (see top panel) the mean nebular intensity distribution cor- 
responding to the averaged flux of the brightest lines H/5 and 
[O III] A5007. The same applies to the continuum emission; no 
significant variation is seen in the intensity distribution from 
4400 A to 6500 A, at least for the small portion of galaxy shown 



sity of the blue and red W-R bumps is observed in nearly all 
the galaxies. This difference is however too small (< 1") to be 
significant, given the spatial resolution of our observations. We 
thus plotted the mean intensity distribution of W-R emission 
lines by averaging the intensities of He II A4686, C IV A5808 
and N III A4640H-C III A4650 (for NGC 3049 and Tol 89 only). 

In the middle panel of each figure, we include an isophotal 
map of the galaxy produced from a CCD image taken with- 
out filter just before the spectroscopic observation. The posi- 
tion and direction of the slit are indicated. This image is very 
useful to identify the peak intensities of the continuum and 
emission lines. In the bottom panel, we show a velocity curve 
derived from the measured central wavelength of nebular emis- 
sion lines along the slit. The mean radial heliocentric veloc- 
ity has been derived using the brightest nebular lines, i.e. H/3, 
[Oiii] A4959, [Oiii] A5007 and Hei A5876, since no system- 
atic variation in the velocity curves has been observed for the 
various emission lines. 

One of the most striking features which appear in Figs. [l| 
and H is the difference between the distributions of stellar con- 
tinuum and nebular lines in He 2-10 and Tol 89. The stellar 
continuum shows only one maximum whereas the distribution 
of nebular lines intensity is mostly double-peaked. Moreover, 
the maximum intensity of nebular lines is shifted relative to the 
maximum intensity of stellar continuum, both for He 2-10 (~ 
1" to the East) and Tol 89 (~ 2" to the South- West). In the other 
galaxies, i.e. NGC 3049, NGC 3125 and NGC 5253, the peaks 
of the nebular lines and stellar continuum intensity distribution 
are coincident. 

In all the galaxies, the peak intensity of the broad emission 
lines He ll A4686 and C IV A5808 coincide with the peak inten- 
sity of the stellar continuum, even in He 2-10 and Tol 89 where 
there is a shift between the distributions of stellar continuum 
and nebular lines. These results confirm the stellar origin of 
the broad W-R emission lines. Given the good signal-to-noise 
ratio of our observations, it has been possible to detect the W- 
R lines over a large extension around the maximum intensity. 
The size of the "W-R regions" ranges from 2" - 3" for Tol 89, 
NGC 3049, NGC 3125 and NGC 5253 to 6" for He 2-10. The 
W-R lines are detected in two distinct regions in NGC 3125 
(separation - 10") and NGC 5253 (separation ~ 3"). See Ta- 
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Fig. 1. He 2-10. Top: Spatial distribution of the stellar contin- 
uum (solid line), of the brightest nebular emission lines (H/3, 
[Olll] A5007, dotted line) and of the broad emission lines 
(He II A4686 and C IV A5808, dashed line) due to W-R stars. 
Capital letters, if present, indicate the position of different WR 
regions (see Table Middle: Isophotal contour map of the 
galaxy showing the position of the slit used for the spectro- 
scopic observations. Bottom: Velocity curve of the ionized gas 
along the slit. The position of peak intensity from W-R stars 
is indicated by the vertical line. Position along the slit is given 
in arcsec. The scale corresponding to 100 pc is given on the 
bottom right corner of the figure 



3.2. Spatial distribution of Wolf-Rayet stars 

In this Section we shall discuss the spatial distribution of the 
W-R signatures in the individual objects. To this end, we use 
(when available) high resolution images from the literature to 
identify the W-R regions. 
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Fig. 2. NGC 3049. See Fig. |l] for the legend 

most prominent is region A at the center of the galaxy. The 
second starburst region B is located at ~ 8.5" East of A. It 
is smaller and fainter in optical images and is seen at an even 
lower level in Ha. The third region called C (Sauvage et al., 
1997), located ^ 2.5" West of region A, is not prominent in 
optical images but is more noticeable in Ha and [O III] images 
( Corbin et al., 1993 ). Regions A, B, and C are resolve d into 
several knots on HSTUY images ( Conti & Vacca, 1994 ). The 
presence of three starburst regions in the central part of He 2- 
10 is consistent with the spatial distribution of W-R emission 
lines along the slit (^ 4" long) which shows a maximum in A 
but extends also to the region C, the secondary peak of neb- 
ular emission seen in Fig.|l]. Knot B was not included in our 
slit. Our long-slit spectrum of He 2-10 shows W-R features in 
regions A and C. However, the spectral signatures in region C 
are too faint to be considered separately from region A in the 
following analysis (see Sect. m. 



He 2-10 



NGC 3049, NGC 3125 



Three starburst regions can be distinguished on the optical im- 
ages of He 2-10 (Hutsemekers & Surdej, 1984; Corbin et al.. 



In NGC 3049 and NGC 3125, the W-R regions coincide with 
the peaks of the continuum emission and nebular lines dis- 
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Fig. 3. NGC 3125. See Fig.[l]forthe legend 



Fig. 4. NGC 5253. See Fig. |l] for the legend 



the Virgo cluster This low-mass galaxy has a small bulge and 
a thin bar of constant surface b rightness surrounded by an 
inner ring ( Contini et al., 1997 ). The optical appearance of 
NGC 3125 is an amorphous elliptical shape with a bright cen- 
tral starburst region dominated by two bright knots (regions A 
and B) apparently connected by a bridge of fainter intensity. 

NGC 5253 

Our observations of NGC 5253 have already been presented in 
Schaerer et al. (1997). The spatial distribution of W-R stars in 
this object is complex (see Fig. There are two W-R regions, 
labeled as A and B, located at the two peaks intensity of the 
continuum and nebular lines distributions. The brightest W-R 
region (B), however, corresponds to the faintest continuum and 
emission lines region. For this galaxy we made a distinction be- 
tween the distribution of the He II A4686 and C IV A5808 lines. 
There is indeed a third bright "He II A4686 region" (region H 
in Schaerer et al., 1997) located 6" North-East of the brightest 
W-R region B. The spectrum of this region show a narrower 
He II A4686 line than in others regions A and B, probably of 
nebular origin (Schaerer et al., 1997). In fact, the core of the 



well as diffusely distributed massive stars (Meurer et al. 1995 
;cf. also Calzetti et al., 1997 for a recent study of NGC 5253). 

Tol 89 

Tol 89 is a giant Hll region located in one of the spiral arms of 
the barred spiral galaxy NGC 5398, about 34" South- West of 
the nucleus. This object shows the largest shift between the spa- 
tial distribution of the continuum and nebular lines. The spatial 
distribution of the W-R lines is narrower than that of the neb- 
ular lines. W-R lines are found at the two regions where the 
continuum peaks. 

4. The Wolf-Rayet regions 

In this section, we present a spectroscopic analysis of the indi- 
vidual regions where spectral signatures from W-R stars have 
been detected. One-dimensional spectra were extracted corre- 
sponding to the different regions by adding several columns 
along the slit. The size of the extraction window has been cho- 
sen for each region in order to obtain the best compromise be- 
tween a high signal-to-noise ratio in the continuum (around 
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Fig. 5. Tol 89. See Fig. [l] for the legend 



We extracted one-dimensional spectra centered on the maxi- 
mum intensity of W-R emission lines with an aperture of ^ 
1.7" for NGC 5253 (regions A and B) and Tol 89, - 2.0" for 
NGC 3049, - 2.7" for He 2-10 and NGC 3125 (regions A and 
B). We thus obtained spectra for a total of seven regions, one 
per galaxy except for NGC 3125 and NGC 5253 in which two 
W-R regions have been identified. The full wavelength range 
(4400 - 6500 A, corrected for the Doppler effect) of the flux 
calibrated spectra (not corrected for reddening) of W-R regions 
are shown in Fig. |[ Normalized spectra showing the region 
around the blue and red W-R bumps are given in Figs. ^ and ^ 
respectively. 

Given the complex structure and multiple blends around the 
blue W-R bump and the medium resolution of our spectra we 
adopt the following procedure to measure accurately the indi- 
vidual emission lines. We use a procedure of multi-gaussians 
fitting where the number of gaussians is the only fixed parame- 
ter (central wavelength, intensity and FWHM are free to vary). 
The normalized spectra were obtained by dividing the spectra 
by a low-order polynomial fit on the line-free region of the con- 
tinuum. In Table we report the measurements performed on 




4500 6000 5600 6000 6500 

Wavelength [A] 

Fig. 6. Optical spectra of W-R regions at the rest frame wave- 
length of galaxies. Note the blue continuum and the bright neb- 
ular emission lines (H/3, [O III], and He I A4471) typical of star- 
burst regions. Normalized spectra showing the region around 
the blue (around 4700 A) and red (around 5800 A) W-R bumps 
are given in Figs. ^ and |^ respectively. The spectra are not cor- 
rected for extinction 



4.7. Uncertainties 

There are three major sources of uncertainties in the measured 
fluxes and equivalent widths of emission lines. The first one 
is inherent to the measurement itself and depends mainly on 
the signal-to-noise ratio of the spectra. These errors were com- 
puted directly when using the IRAF task SPLOT to derive 
the parameters of emission lines. The second source of errors 
comes from the extraction of an individual spectrum corre- 
sponding to a region in a galaxy. Indeed, for some galaxies like 
He 2-10 and Tol 89, we observe a shift between the maximum 
emission of the stellar continuum and of the nebular lines (see 
Sect. 1^). Thus, the measured line flux and equivalent width of 
the lines is sensitive to the size of the extraction window used 
to obtain one-dimensional spectra of the individual regions. To 
estimate this error, we measured the lines in extracted spectra 
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Fig. 7. Enlargement of the blue W-R bump in the optical spec- 
trum of galaxies. Spectra are normalized to the continuum and 
an offset of 0.4 has been applied between each spectrum. The 
broad He II A4686 line from WNL stars is clearly detected. The 
broad line around 4645 A is not well identified and could be a 
blend of N III A4640 and C III A4650 lines due to WN and WC 
stars. The spectra of NGC 3 125 -A and NGC 5253- A exhibit the 
typical high-excitation nebular emission lines of [Ar IV] A47 1 1 
and [Ar IV] A4740. Their presence might imply a nebular con- 
tribution to the He II A4686 line 

source of uncertainty exists for lines from the blue W-R bump. 
Given the important slope of spectra and the complex struc- 
ture in this region the line measurement also depends on the 
determination of the continuum level (i.e. the normalization). 
Therefore measurements of emission lines in the W-R bumps 
were performed before and after normalization. 

The three sources of errors were propagated quadratically 
to obtain the final uncertainties reported in Table || for the 
dereddened line fluxes and for the equivalent widths. For the 
bright nebular emission lines (like H/3, [O III] and He I A5876), 
the relative uncertainty is generally smaller than 10%. How- 
ever, it can reach a maximum of ~ 20% for the equivalent 
width and of ~ 30% for the relative flux of faintest broad emis- 




5600 5700 5800 5900 

Wavelength (A) 



Fig. 8. Enlargement of the red W-R bump in the optical spec- 
trum of galaxies. Spectra are normalized to the continuum and 
an offset of 0.2 has been applied between each spectrum. The 
broad C IV A5808 line from WC stars is clearly detected in all 
spectra. Except for He 2- 1 0, this represents the first detection of 
WC stars in these galaxies. Note also that a broad emission line 
of C III A5696 is also reported at ^ 2a level in the spectrum of 
NGC 3049. This indicates the predominance of late-type WC 
stars in this high-metallicity galaxy 



the measurement itself (which is cx S/N) and to the extraction 
window are comparable. For emission lines from the blue W-R 
bump, the uncertainties due to the normalization are not negli- 
gible and can even be comparable in some cases. We note how- 
ever some exceptions. For He 2-10 and Tol 89, which show a 
discrepancy between the nebular and the continuum emission, 
the uncertainty on the equivalent width of nebular lines (and 
especially for H/3) is dominated by the error due to the size of 
the extraction window. 

We checked the absolute flux calibration (H/3 flux), the rel- 
ative line flux of He I A4471, [Olll] A4959 and He I A5876 
lines, and the equivalent width of H/3 by comparing our mea- 
surements with previous observations reported in the literature: 
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1986; IVacca & Conti, 1992|; Pontini et al., 1997| ), NGC 3125 
dKunth & Sargent, 198 1|; ^Cunth & Sargent, 1983| ; Kunth & Jou- 
bert, 1 985; |Vacca & Conti 1992|), NGC 5253 (Ca mpbell et al., 
1986; [Walsh & Roy 1989|; tTerlevich et al., 199 1| ), and Tol 89 
( Durret et al., 1985 ; rerlevich et al., 1991 ). A good agreement 
was generally found for the absolute calibration and the relative 
fluxes of bright nebular lines, with relative errors ^ 20%. 

Larger differences can be found for the broad W-R lines. 
Comparing our data with VC92, our measurements of Hell 
A4686 (relative flux to H/3 and equivalent width) are found to 
be systematically larger by a factor of ~ 3 in NGC 3049, and 
~ 1.5 in He 2-10 and NGC 3125, although a good agreement 
is found for the nebular lines. Given the systematic differences, 
it is unlikely that different slit positions are responsible for this 
effect. To understand the origin of this discrepancy, we com- 
pared our spectra to those published in VC92 and kindly pro- 
vided by William Vacca. The comparison shows that the differ- 
ences can mostly be understood in terms of differing signal-to- 
noise ratios. Indeed, the measurement of the faint, broad W-R 
emission lines is strongly dependent on the signal-to-noise ra- 
tio, especially to fit the continuum level in the deblending pro- 
cess used to separate the W-R bump into the individual compo- 
nents. Given the improved quality of our spectra (S/N ^ 50-70 
around 4700 A compared to S/N < 40 for VC92), our mea- 
surements should be more reliable. In particular the improved 
S/N also allows the detection of the weaker N III A4640H-C III 
A4650 blend, barely detectable in the data of VC92. 

Within the uncertainties, our new measurements of W-R 
emission lines in NGC 5253 are consistent with those given 
in Schaerer et al. ( 1997 ), which were obtained from a different 
reduction/analysis. 

4.2. Reddening 

Due to the limited spectral range, our observations do not in- 
clude the Ha emission line commonly used to determine the 
internal reddening parameter ch/3. We therefore adopt values 
taken from the recent literature. We wish to remind that the de- 
termination of the amount of reddening is subject to two com- 
plications often not accounted for: 1) the value of the extinc- 
tion coefficient can be wavelength dependent, and 2) variations 
of the extinction on small spatial scales can occur. The former 
is well illustrated by the case of He 2-10, where ch/3 varies 
from ^ 0.4 in the optical (Sugai & Taniguchi, 1992; Vacca & 
Conti, 1992) to - 4.5 in the infrared (lAitken & Roche, 19841; 



Kawara et al., 1989 ) and can even reach ~ 14 in the millimeter 



range (Kobulnicky et al., 1995). NGC 5253 is a good example 
showing large spatial variations of reddening (Walsh & Roy, 
1989; Calzetti et al., 1997). The adopted extinction coefficient 
is taken from regions corresponding as closely as possible to 
the ones studied here and derived from optical spectra, i.e. from 
the Balmer decrement. Our adopted reddening values are given 
in Table |[ 

For He 2-10 and NGC 3125, we adopt the internal extinc- 
tion coefficients given by VC92, with different values for re- 



NGC 3049, we cho ose to use the value ch/3 ^ 0.23 given by 
Contini et al. ( |l997| ) instead of that derived by VC92. Our W-R 
region is indeed better identified in the spectroscopic observa- 
tions of Contini et al. ( 1997 ). For Tol 89, we adopt the redden- 
ing coefficient ch/3 ~ 0.18 estimated by Terlevich et al. ( 1991 ) 
which is in good agreement wit h the m easurements of Durret 
et al. ( |1985| ). See Schaerer et al. ( |1997| ) for NGC 5253. 

Dereddened fluxes are derived using the extinction law of 



Cardelli et al. (1989) and including a Galactic foreground ex- 
tinction E{B — V) given in Table [l|. The result is almost in- 
dependent of the chosen extinction curve, since in the wave- 
length range 2600 — 9000 A, there is little difference among 
the mean interstellar extinctions of the Galaxy, the Large Mag- 



ellanic Cloud, and the Small Magellanic Cloud (Seaton, 1979 



Bouchetetal., 1985 ; Fitzpatrick, 1986 ). Dereddened line fluxes 
are given in Table y. 



4.3. Oxygen abundance 

Oxygen is an important element for our subsequent analysis, 
as it is used to define the metallicity of each W-R region for 
the comparison with starburst models (see Sect. ||). Since the 
[O II] A3727 and [O III] A4363 lines are not included in our 
wavelength range, the oxygen abundance was determined from 
t he empirical calibration between log(0/H) and the ratio i?3 
( [Edmunds & Pagel, 1984|) , with the linear- fit given by VC92. 
Edmunds & Pagel (1984) estimate the intrinsic uncertainty 
in the calibration to be about ±0.2 in log(0/H). For the so- 
lar oxygen abundance, we adopted the local galactic value of 
log(0/H)Q = -3.08 ( [Meyer, 1985[ ). The oxygen abundances 
derived in this way are summarized in Table ^. Except NGC 
3049 which has an oxygen abundance slightly above solar, the 
values for the other W-R regions are below solar, typically 
0.2- 0.5 Zq. 

Our derived oxygen abundances are generally in good 
agreement with earlier studies. The oxygen abundance (O/H) 
^ 0.2 Zq derived for NGC 5253 is identical to the value ob- 
tained by Kobulnicky et al. (|1997[). We obtain (O/H) - 1.2 Zq 



for the W-R region of NGC 3049, in good agreement with pre- 
vious observations by Contini et al. (1997). A good agreement 
is also found for NGC 3125 between our measurements (0.17 
Zq for region A and 0.20 Zq for region B) and those of VC92. 
The same is true for Tol 89, when we compare our estimation 
of (O/H) - 0.25 Zq to the derivation by Durret et al. ( |1985[ ). 

For He 2-10, however, we obtain (O/H) ^ 0.42 Zq, three 
times larger than the value from VC92 for the same region. 
The calculations by VC92 are based on the detection of the 
temperature-sensitive [Olll] A4363 emission line, which is 
however very faint in their spectrum of He 2-10. Their (O/H) 
value may therefore be underestimated, and we thus prefer the 
higher value derived from our spectra. 

4.4. Broad emission lines from Wolf-Rayet stars 

Several broad emission Unes (FWHM ^ 25 A) due to W-R 
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Line/Galaxy 


He 2-10 


NGC 3049 


NGC 3125 


NGC 5253 


Tol 89 








A 


B 


A 


B 




Nebular Emission Lines 


TJa T \ AAl 1 

rlc i / i 


1 1 -1- "5 
i lit J 


1 1 -1- 


^Zzt o 




AQ4- 1 


4 J It J 


^ 1 -1- 1 Q 


[re llij y\'4-0JO 


JJ± J 












97-1- S 
Z /zt o 










1 i It i 


1 A4- A 






TArTVl A4740 






4± 1 


3± 1 


8± 2 






He I A4922 






7± 2 


8± 2 


9± 1 


9± 1 


11± 2 


[Oiii] A4959 


471 ±22 


106±6 


1970± 67 


1526± 41 


2070± 36 


1491±19 


1229± 33 


[Oiii] A5007 


1430±44 


338±4 


5852±197 


4538±128 


6155±104 


4389±59 


3636±111 


[Ni] A5199 


12± 2 


14± 5 




7± 2 


9± 1 


8± 2 


3± 1 


[Felll] A5271 










5± 1 


8± 1 




[CI III] A5518 






6± 1 




4± 1 


5± 1 




[CI III] A5538 










3± 1 






[Nil] A5755 


10± 2 


8± 1 






5± 1 






He I A5876 


123± 7 


133± 2 


121± 3 


120± 6 


131± 2 


130± 1 


118± 3 


[Oi] A6300 


14± 5 


23± 1 


18±3 


19± 1 


27± 5 


28± 2 


13± 1 


[Siii] A6318 


21± 7 


6± 1 


14± 1 


12± 1 


29± 7 


20± 1 


20± 1 


[0 1] A6364 


4± 1 


7± 1 


5± 1 


7± 1 


7± 1 


8± 1 


7± 4 


am 


196± 8 


16± 3 


66± 11 


73± 8 


212± 7 


49± 2 


13± 1 


W{R(3) [A] 


23± 3 


34±3 


93±9 


70±7 


269± 23 


112± 10 


68± 17 


Ch/3 


0.56 


0.23 


0.40 


0.64 


0.44 


0.20 


0.18 


(O/H) [Zg] 


0.42 


1.20 


0.17 


0.20 


0.20 


0.20 


0.25 


Broad Emission Lines from Wolf-Rayet Stars 


N III A4640+C III A4650 


112±43 


117±35 


31±4 


63± 18 


7± 2 


18± 8 


81± 32 


Hen A4686 


97±35 


128±38 


76± 8 


73± 15 


9±3 


45±15 


116± 43 


Cm A5696 




35±10 












Civ A5808 


41± 7 


79±10 


21± 5 


51±9 


6±3 


30±5 


132± 30 


W(NIII A4640+C III A4650) [A] 


2.7±0.5 


5.2±0.5 


3.2±0.5 


2.8±1.0 


2.1±1.3 


1.0±0.5 


3.4±0.5 


W(Heii A4686) [A] 


2.3±0.5 


5.6±0.5 


7.3±1.0 


2.8±0.5 


3.0±0.7 


3.3±1.1 


6.1±0.5 


W(Cm A5696) [A] 




2.0±0.5 












W(Civ A5808) [A] 


1.5±0.5 


4.7±0.6 


3.1±0.4 


6.2±1.0 


2.5±0.6 


6.0±1.1 


11.9±1.5 



For each emission-line we reported the dereddened flux normalized to I(Hl 3) = 1000. /(H/9) is the absolute dereddened flux of the H/3 



emission line in 10 ergs s cm , ch/3 is the extinction coefficient (see Sect. 4.2 for references) and W\ is the equivalent width of the H/3 



and W-R emission lines in A. (O/H) is the oxygen abundance in solar units (see Sect. 4.3 for details). 



over the wavelength range covered by our spectra: the blend 
of N III A4640 and C ill A4650, He ii A4686, and C IV A5808. 
Equivalent widths and fluxes relative to H/3 of these emission 
lines are listed in Table |[ 

One of the main results of this paper is the unambiguous 
detection of broad (FWHM - 50 - 90 A) C IV A5808 emission 
in all W-R regions of the observed galaxies (see Figs. |] and 

which clearly indicates the presence of WC stars in these 
regions (cf. below). Whereas all the objects where previously 
known to show the broad W-R bump around 4700 A, our spec- 
tra are the first to show the red W-R bu mp aro und 5800 A. 
As already pointed out by Schaerer et al. ( 1997) a sufficiently 
high signal-to-noise ratio (> 40) is required for the detection 
of such a C IV A5808 line in integrated spectra of W-R galax- 
ies. Typically the 5808 flux is found to be factor of two less than 
4686; furthermore the larger width of C IV A5808 and the lower 
intensity of the continuum render its detection more difficult. 



ous ground-based spectroscopical studies, except for He 2-10 



(Hutsemekers & Surdej, 1984; Vacca & Conti, 1992) 



What do the broad emission features tell us about the W-R 
populations in these regions ? WN stars cannot be responsible 
for the C IV A5808 emission, since they show He II A4686/C IV 
A5808 ~ 16. (SV98), much larger than observed in our spec- 
tra. In addition the strong He II A4686 emission, the rela- 
tive strengths of the W-R features in the blue bump, and the 
FWHM(He II A4686) (- 30 - 40 A) are fin gerprints of the 
presence of WN stars, as commonly accepted ( |Conti, 1991 ). In 
all our seven emission line regions we thus find signatures of 
both WN and WC stars. For six regions this represents the first 
detection of WC stars. 

The dominant W-R subtypes can be constrained as follows. 
In all our spectra the presence of N III A4640 and/or C III A4650 
is established, while N V A4604 is very weak or absent. If the 
former are due only to WN stars this indicates the predomi- 
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(Nlll A4640+CIII A4650)/Hell A4686 ranges from 0.9 to 2.5 
with an average of 1.5, which is much larger than the average 
of 0.24 ± 0.15 obtained for LMC WNL stars (SV98). The ex- 
cess flux can easily be attributed to the emission from the WC 
stars. As for the WN stars it is not possible to apply the usual 
criteria to classify the WC subtypes. However, the measured 
FWHM(5808) are comp atible with types W C7-WC4 but also 
WOl, W03, and W04 ( |Smith et al., 1990| ; Crowther et al., 
1998). A significant population of WOl stars can, however, be 
excluded from the absence of O V A5590. Only in NGC 3049 
do we also detect C III A5696 (cf. below) with a relative in- 
tensity compatible with WC7 stars. In the other cases the non- 
detection of C III A5696 and FWHM(5808) favors early-type 
WC stars. From the comparisons of the broad emission fea- 
tures, we thus conclude that all regions (except the high metal- 
licity galaxy NGC 3049 discussed below) contain a mixed pop- 
ulation of WNL, and early WC and/or W03-4 stars. The pres- 
ence of intermediate-type WNAVC stars cannot be excluded. 

In the following we shall now discuss the W-R signatures 
individually for each object. 

He 2-10 

He 2-10 w as the first galaxy discovered to have evidence for 
W-R stars ( Allen et al., 1976 ), a nd is in many ways considered 
as the prototypical W-R galaxy ( Conti, 1991 ). A broad feature 
of C IV A5808 was also suspected by Hutsemekers & Surdej 
(1984) and detected by VC92 in the optical spectrum of He 2- 
10. The presence of this line is unambiguous in our spectrum 
of He 2-10 (see Fig. The blue W-R bump (see Fig. ^ is 
composed by the broad He II A4686 fine and N III A4640H-C III 
A4650 with nearly the same strength. 



NGC 3049 

In NG C 3049, the W-R region first noted by Kunth & Schild 
( 1986 ) is located in the bar, ^ 2.5" North-East of the galaxy nu- 
cleus (see Fig. H). The blue W-R bump is defined by two broad 
emission lines identified as He II A4686 and N III A4640H- C III 
A4650 (see Fig. Q), with approximatively the same strength 
as already reported in previous observations (Kunth & Schild 
1986, VC92). In addition to the CiV A5808 fine, a broad 
(FWHM 45 A) emission line of C III A5696 is also detected 
at 2a level in the specti 'um of NGC 3049 (see Fig. |). Af- 
ter Phillips & Conti (1992) this is the second detection of this 
W-R emission line in an extra-galactic Hll region. The C III 
emission is attributed to late WC stars, whose presence is in- 
deed expected in high metallicity regions (Maeder, 1991 ). 



NGC 3125 



Penston et al. (1977) detected the He II A4686 emission line in 
the North- West knot (= our region A) of NGC 3125, but do 
not mention W-R stars. The W-R nature of NGC 3 125 was first 



cesses in Hll galaxies. They computed that this starburst region 
contains several hundred W-R stars. Long-slit observations of 
VC92 have first revealed W-R stars in the second condensation 
10" South-East of the nucleus (= our region B). In region A of 
NGC 3125, the W-R bump around 4700 A is dominated by the 
He II A4686 line whereas the N III A4640H- C III A4650 blend 
and the He II A4686 line are of similar strength in region B (see 
Fig. We signal a marginal la) detection of broad He II 
A5412 emission, the second strongest Hell line in WN stars, 
in knot A. If confirmed this represents, to our knowledge, the 
first detection of this line in the integrated spectrum of an extra- 
galactic object. The spectrum of region A also exhibits the typi- 
cal high-excitation nebular emission lines of [ Ar IV] A47 1 1 and 
[Ar IV] A4740 (see Fig. Their presence might imply a neb- 
ular contribution to the He II A4686 line. 

NGC 5253 



The spectra are discussed in Schaerer et al. ( 1997 ). In addition 
we signal that Hell A5412 emission may also be marginally 
detected in region B. 

Tol 89 



Durret et al. ( |1985| ) reported the first detection of the W-R bump 
around 4700 A in the spectrum of Tol 89, but their low signal- 
to-noise ratio and moderate spectral resolution did not allow a 
clear identification of the lines. They do not detect any optical 
spectral signature from WC stars, but their presence was sug- 
gested from UV P-Cygni lines. In our spectrum of Tol 89 (see 
Fig. the blue W-R bump shows He II A4686 and at a lower 
level the N III A4640H- C III A4650 blend. The C IV A5808 fine 
is the strongest (14^ ~ 13 A) and the broadest (FWHM 90 
A) detected among the spectra of observed W-R regions (see 
Fig.|). 

4.5. Massive star populations 

In this section we derive the approximate number of massive 
stars (O and W-R) in the W-R regions of galaxies. The me- 
chanical energy output from the stellar population and its im- 
pact on the surrounding ISM is studied in Sect. A direct 
comparison of the observed W-R emission line features with 
evolutionary synthesis models is presented in Sect. |[ 

The approximate number of WN, WC, and O stars are listed 
in Table g (columns 4, 5 and 7). The number of W-R stars is 
calculated from the luminosity of the W-R emission lines. We 
assume that the dominant contributors to the broad He II A4686 
and C IV A5808 lines are respectively WNL and WC4 stars (see 
Sect. 4.4). The average observed luminosity of WNL stars in 



noted by Kunth & Sargent (1981 ) who first envisioned thekim 



the He II A4686 Une is L4686 = 1-6 ± 1.5 x 10^*' ergs s"^ fliat 
of WC4 stars in the C IV A5808 line is L58O8 = 3.0±1.1 x lO^e 
ergs s^^(SV98). The quoted uncertainty on the WN num- 
ber reflects the standard deviation of i4686- Given the com- 
plex structure of the blue W-R bump, the uncertainty on He 11 
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Table 3. Massive star populations and mechanical energy input rates in W-R regions 



Galaxy 




Age 
[Myr] 


Qo 
[10*« s-M 


NwN 


iVwc 


Vo 




No 


log AL 

[Mq] 


logs 
[ergs s"^] 


[pc] 


Vb 

[kms-^] 


(1) 




(2) 


(3) 


(4) 


(5) 


(6) 




(7) 


(8) 


(9) 


(10) 


(11) 


He 2-10 




5.5-6.0 


3802 


1100±520 


> 250 


0.5-1.0 


2450 


-4900 


6.8 


40.81 


1000 


110 


NGC 3049 




5.5 


1349 


510±240 


> 170 


a 






6.4 


40.40 


900 


100 


NGC3125 


A 


4.5-5.0 


2188 


500±230 


> 70 


0.25-0.5 


3240 


-6470 


6.1 


40.48 


700 


90 




B 


4.5-5.0 


2455 


530±250 


> 200 


0.25-0.5 


3450 


-6900 


6.1 


40.53 


700 


90 


NGC 5253 


A 


3.0 


896 


26±13 


>9 


0.8-0.9 


960 


- 1080 


6.6 


39.00 


230 


50 




B 


5.0 


207 


27±13 


> 10 


0.25 




680 


6.6 


39.30 


440 


50 


Tol 89 




4.5-5.0 


708 


240±110 


> 150 


0.25-0.5 


640 


- 1270 


5.7 


39.84 


550 


70 



° For instantaneous burst models with solar metallicity and higher, rj is not defined. The derived number of O stars (col. 7) is strongly 
dependent on the adopted temperature limit between O and B stars. See text for more details. 



lar He II contribution may also be present in some of our ob- 
jects (NGC 3125-A, NGC 5253-A, and possibly also Tol 89; 
see Fig. In this case iVwN has to be corrected downwards. 
Since L58O8 is lower for later WC types (SV98) and the dom- 
inant WC type is not well constrained (see Sect. p~4| ) a lower 
limit on -/Vwc is provided. Taking into account the variation in 
-L58O8 between WC4-6 stars we conclude that the number of 
WC stars may be larger by a factor of ~ 3 — 5. 

Under the condition of case B recombination and assum- 
ing that all the ionizing photons emitted by the stars are ab- 
sorbed by the gas, the total number of Lyman photons, Qq^^ 
(col. 3), can be derived from the observed luminosity of the 
H/3 emission line. Note that this also provides a convenient 
estimate of the number of equivalent 07V stars, for which 



Q^''^ = 1.0 X 10*3 s-i dLeitherer, 199^ . To estimate the 



number of O stars, A^o^ we must take into account the ioniz- 
ing photon contribution from W-R stars, the age of the stellar 
population and the IMF (see S V98). iVo is given by: 



Nr 



/^obs 



AT ^W-R 



,07V 



(1) 



where A'^vv-R — Nwn + Nwc^ Q'q ^^ is the average Lyman 
continuum photon flux per W-R star, and rjo {t) is the IMF aver- 
aged ionizing Lyman continuum luminosity of a stellar popula- 
tion normalized to the output of a 07V star (see SV98). In most 
regions we find A^vv-r/^07V = (A^wn + Nwc)/No7y ^ 
0.2 — 0.5. Assuming a priori that the massive main sequence 
stars provide the bulk of the ionizing Lyman continuum pho- 
tons implies that the average Lyman flux per W-R star, Q^^, 
is not significantly larger than Q^^^. For simplicity we adopt 
qW-R ^ Q07V ^ 1.0x10^9 s-i (similar to VC92 but smaller 
than the age dependent value given by SV98). The value of 770 
given in column 6 of Table |] was taken from the instantaneous 
burst models of SV98 for a Salpeter IMF at the age given in 
column 2 (derived from W{Il(3)). The resulting number of O 
stars (defined by Tcs > 33 000 K) is given in column 7. The in- 
dicated range results only from the range of 770. The number of 
O stars found for the W-R regions are between ~ 500 and 7000. 



age (determined from W{l[lf3)) the instantaneous burst models 
at solar or higher metallicity predict that no stars with Tcff > 
33 000 K are left. If this temperature limit is adopted to separate 
O stars from later types (cf. S V98) we formally derive A^o = 0. 
(see Table ||). Lowering this limit by 10 % already changes the 
situation qualitatively: in this case one obtains rj ^ 0.16, and 
A^o ^ 4000 as test calculations show. In any case, according to 
these models, the bulk of the ionisation is provided by late O 
and/or B type stars and W-R stars. 

From Table |[ we obtain W-R/O number ratios of ^ 0.03 — 
0.6 (and even larger for NGC 3049), systematically higher than 
the predictions for constant star formation at the appropriate 
metallicity (Maeder & Meynet, 1994), but within the range of 
instantaneous burst models with different IMF slopes (SV98, 
Table 4). A more detailed attempt to constrain the IMF is made 
in Sect. 5.3. A trend of increasing W-R/O ratios towards higher 
metallicity is found as expected on the average (e.g. Meynet 
1995). 

For the seven W-R regions, we find a flux ratio for He II 
A4686/C IV A5808 of 1 .8 ±0.8, which shows a fairly smafl dis- 
persion. The corresponding number ratio of WC/WN stars is 
typically 0.2 - 0.4; the extreme values are 0.14 (NGC 3125- 
A) and 0.63 (Tol 89). No systematic variation with metallicity 
is found. For the low-metallicity objects (Z ^ 0.2 Zq; NGC 
3125, NGC 5253, and Tol 89) the derived WC/WN ratio is 
larger than what is found in the Local Group (except IC 10) at 
similar metallicity ( ^assey & Johnson, 1998 ). At higher metal- 
licities He 2-10 and NGC 3049 show, however, WC/WN ratios 



below the trend observed by Massey & Johnson ( 1998 ) 



The finding of a fairly constant WC/WN ratio may seem 
surprising at first sight. Contrary to regions of constant star 
formation, likely representative of the Local Group samples, 
regions of short star formation (more appropriate for our ob- 
served starburst galaxies, see Sect, ^j] ) could a priori show 
quite a large range in WC/WN, depending on the age of 
the starburst. This is illustrated in Fig. |^ where the predicted 
WC/WN ratio from the SV98 models are shown for different 
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Meynet, 1995[])- However, already burst durations of ^ 2 Myr 
suffice to smooth out the rapid variations of WCAVN. This may 
well explain the small range of WC/WN found for the low- 
metallicity objects. The observed WC/WN value in these ob- 
jects is also intermediate between the predictions from high and 
standard mass loss models (cf. Fig. ^a and^). More surprising 
is the low WC/WN ratio of the remaining higher metallicity ob- 
jects (He 2-10 and NG C 3049 ) compared to the WC/WN value 
of Massey & Johnson (1998) at a similar metallicity. Indeed, 
the probability of finding WC/WN ratios below the equilibrium 
value attained in regions of constant star formation should be 
quite small. However, if the number of WC stars is system- 
atically underestimated by a factor of 3 (due to variations of 
the average WC subtype with metallicity; cf. above) our ob- 
servations may well all be larger than the observed WC/WN 
trend with Z. We conclude that the WC/WN ratios of the low- 
metalHcity galaxies NGC 3125, NGC 5253 and Tol 89 can 
be understood quantitatively with burst models of reasonably 
short but non-zero duration. Additional observations of WC 
and WN populations, especially for regions of higher metal- 
licities, would be very helpful. 



4.6. Energy released by massive stars 

The winds produced by these populations of massive stars and 
by supernovae explosions should strongly perturb the gas dy- 
namics of galaxies. We effectively observe some perturbations 
centered on the W-R regions in the velocity curves of the ion- 
ized gas (see Figs. |l] to |^. The observed velocity differences 
range from - 40 km s^^ in NGC 3125 and NGC 5253 to - 
100 km s^^ in He 2-10 and NGC 3049. These values are similar 
to the typical expanding velocities (~ 50 km s^^) of superbub- 
bles and filaments observed in dwarf galaxies (Marlowe et al., 
1995; Legrand et al., 1997a). According to these authors, the 
mechanical energy output from the supernovae explosions and 
strong stellar winds in the starburst regions appears adequate 
to power expansion motions of this speed. Note also that gas 
flows resulting from the energizing of the interstellar medium 
are currently observed in star-forming galaxies (Kunth et al., 
1998). 



To verify this claim, one can compare the computed expan- 
sion velocity of a bubble of gas centered on the W-R regions to 
the observed velocity differences in galaxies. We used a very 
simple model of wind-blown bubble expanding adiabatically 
(i.e. e nergy conservative), which has been described by Castor 
et al. ( 1975 ) who give the following expressions for the evolu- 



As pointed out in SV98 due to a different interpolation technique 
the models of Meynet ( 1995) predict a more important WC population 
at later ages than SV98 even if the same tracks are used. The "real" 
number of WC stars is likely intermediate between the two models. 
For a given set of tracks the WC/WN plotted in Fig. p| may thus be 
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Fig. 9. Predicted WC/WN number ratio as a function of age 
for stellar evolutionary tracks assuming high mass loss rates 
(a) and standard mass loss rates (b) as described in Meynet 
et al. (1994). Instantaneous burst models for solar metallic- 
ity (Z=0.020, thin dashed Une) and Z=0.008 (thin dotted) are 
shown in a. Thick lines show models at Z=0.004 for differ- 
ent burst durations: instantaneous burst (long-dashed). At = 2 
Myr (dashed-dotted). At — A Myr (dashed), and constant star 
formation (solid). All models assume a Salpeter IMF. 



tion of the bubble radius (in kpc) and expansion speed (in km 

s-i): 




rh 



or, / ^41 \ ,-2/5 
\ "0 / 



(2) 



(3) 



where £^41 = {dE/dt)ii is the kinetic injection rate in units of 
lO''^ ergs s^^, no is the number density (in cm^"^) in the ambi- 
ent interstellar medium, and is the time since the expansion 
of the bubble began in units of 10^ yr. The rate of kinetic energy 
produced by the population of massive stars En can be esti- 
mated from starburst model predictions (Leitherer & Heckman, 
1995) assuming an instantaneous burst and a Salpeter IMF. The 
predicted normalized rate of mechanical energy at a given age 
is thus multiplied by the total mass of ionizing stars formed in 
the starburst region. The total burst mass (column 8 of Table ||) 
is derived from Q'^^ using the instantaneous burst models of 
SV98 and assumi ng a Salpeter IMF dow n to 0.8 Mq. We adopt 
no ~ 0.3 cm~'^ (Marlowe et al., 1995), and for tj, we use the 
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computed radius rb and expansion velocity Wb of the bubble are 
given in the last two columns of Table ^. 

Comparing these values to the velocity differences ob- 
served in the W-R regions of galaxies (from ^ 40 to 100 km 
s^^), we find that the observed and predicted bubble expan- 
sion speeds are in satisfactory agreement, especially when con- 
sidering the extreme simplicity of the model and the observa- 
tional uncertainties. The comparison of radius is more doubtful 
since, with long-slit spectra, we have some informations in only 
one spatial direction. The rough assessments of "bubble" radius 
from velocity curves are systematically lower by a factor two 
(NGC 3049) to six (He 2-10) than the values predicted by the 
model. The sa me trend has been obse rved using Ha images of 
dwarf galaxies ( Marlowe et al., 1995 ). These discrepancies are 
certainly due to the difficulty to measure a "radius" from the ve- 
locity curves and to the simplicity of the model which assumes 
a pure spherical bubble expanding adiabatically. When look- 
ing deeply in Ha images of nearby galaxies, like NGC 5253 
( [Calzetti et al., 1997 ), one can see that the structure of expand- 
ing gas is much more chaotic, ranging from large scale (~ 1 
kpc) filaments to roughly circular shell nebulae of smaller size 
(< 100 pc). It is thus very difficult and even meaningless to 
compare some predicted and observed bubble "radius". 

Nevertheless, the size (~ 50 to 300 pc) of expanding gas 
structures observed in W-R regions is significantly smaller than 
those derived from Ha images (^ 700 to 1500 pc) of dwarf 
galaxies (Marlowe et al., 1995). This might reflect the different 
time scales of the gas dynamics, which depend on the age of 
the starburst responsible for the release of kinematical energy 
through stellar winds and supernovae explosions. In the W-R 
regions, the starburst is one order of magnitude younger (a few 
10^ yr) than the one (a few 10^ yr) derived from gas kinemat- 



ics in older starburst regions ( Marlowe et al., 1995 ). This could 
explain the smaller sizes of expanding "bubbles" in W-R re- 
gions. An additional effect may be the smaller energy injection 
expected from young regions where the SN rate is still small 
(Leitherer & Heckman, 1995). In any case it is evident that 



the large population of massive stars derived in W-R regions 
strongly perturb the surrounding ISM through the injection of 
their mechanical energy. 

5. Comparison with starburst models 

All the observed W-R features and the most important nebu- 
lar H and He lines have been modeled recently with synthesis 
models by SV98. The aim of this Section is to constrain the 
main burst parameters (age, IMF, star formation history) as far 
as possible from the observed W-R signatures. Finally the suc- 
cess or failure to reproduce the observed features also provides 
a test of the underlying evolutionary tracks. 

5.1. Model parameters and comparison procedure 

In the following we consider two basic free model parame- 
ters: the IMF slope and the duration of the star-forming event. 



the observations (see values given below). A power law with 
a slope a is adopted for the IMF (in our notation a = 2.35 
for a Salpeter IMF). The upper mass cut-off is generally set 
to Mup = 120 Mq; the results discussed in this work are not 
affected by the choice of the lower mass cut-off. We consider 
burst models with different durations At starting at time t = 0, 
with the limiting case of an instantaneous burst (At = 0). 

For the subsequent detailed comparison it is, however, use- 
ful to remind some assumptions or uncertainties of the models, 
which are not considered as proper free parameters. 

(a) Most important is the set of stellar evolution models 
adopted. We use SV98 models based on the latest Geneva 
stellar evolution tracks (high mass loss models of Meynet 
et al., 1994), which have been extensively compared to ob- 
servations (see Maeder & Meynet 1994) and which in par- 
ticular reproduce well massive star populations in the Lo- 
cal Group. Only single star models are considered (but see 
Sect. 5.3 for binary stars). 

(b) The predicted fraction of WC stars (WC/WR, and 
WC/WN) is affected by the choice of interpolation tech- 
niques as briefly discussed in Sect. 4.5 and SV98. During 
the WR phase at ages t > 4 Myr the relative WC popula- 
tions predicted by the Meynet (1995) and SV98 models dif- 
fer; compared to the models used here (SV98), more WC 
stars are likely to exist att ^4 Myr. 

(c) As mentioned in Sect. 4.5 the predicted strength of the W-R 
lines is affected by the uncertainty of the intrinsic W-R line 
luminosity, which is the largest for He II A4686. 

(d) To relate models with observations M^(H/3) is conveniently 
used as an age indicator. This assumes that the models pre- 
dict the correct number of ionizing photons (i.e. correct at- 
mospheres, tracks). Case B (ionization bounded nebula) is 
valid and the fraction of photons absorbed by the gas is 
correct (SV98 adopt = 1), and the continuum light is 
correctly predicted. All except the last point determine ob- 
viously the correctness of the nebular line intensity, whose 
strength is also used for comparisons with the WR lines 
(see below). It is fairly difficult to verify how well the pre- 
dicted H/(H/3)-age relation (cf. Copetti et al.l986) holds. 
Empirical tests using recent models and observations of 
Hll regions and their stellar content in the Galaxy and the 
Magellanic Clouds could be useful. The fact that very few 
Hll regions with observed Vl^(H/3) as large as predicted by 
synthesis models for ages ~ 0-2 Myr may indicate a diffi- 
culty in the W^(H/3)-age relation (e.g. Mas-Hesse & Kunth 
1998). 

The comparison with the observations can be performed in 
two ways. One may compare relative line intensities of the W- 
R features with respect to the nebular H/3 emission (in short 
W-R/H/3), or use equivalent widths of the W-R lines. Obser- 
vationaly the determination of these quantities can mainly be 
"perturbed" by three effects: 

( i) We do not count all the ionizing photons produced by mas- 
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Hll region of interest is not entirely included in the slit (e.g. 
if the gaseous component is more extended than the ion- 
izing star clusters), and/or if ionizing photons can escape 
from the Hll region (i.e. not ionization bounded nebula). 

( ii) Stars and gas suffer from a different extinction. 

( Hi) An underlying older population contributes additional con- 
tinuum light. 

The relative intensities W-R/H/3 are effected only by (i) and 
(ii). Effect (i) increase the relative WR/H/3 intensities. The 
same holds for ( ii) if the stellar light is less extincted than the 
gas (cf. e.g. Calzetti 1997, Mas-Hesse & Kunth 1998). The 
equivalent widths are only affected by (Hi), which decreases 
the observed value. For the present paper the quanUtative im- 
portance of these effects cannot be asserted in general without 
additional observations. At least our spatial information (see 
Sect. 3^ ) can provide hints on ( i). From the work of Mas-Hesse 
& Kunth (1998) we see that, if present, the differences in ex- 
tinction between the gas and stars are typically AE{B — V)^ 
0.05 - 0.1. This implies that the stellar flux may be overesti- 
mated by a factor of - 1.2 - 1.4, i.e. /(W - R)//(H/?), iVwN 
and hence -/Vwn / O have to be reduced by the same amount. 
Although e.g. the WCAVN ratio is in principle also modified, 
the effect is small (-4-8 %). For NGC 3 125, NGC 5253, and 
NGC 3049 an underlying population contributes to — 60 % of 
the light at 4630 A (Mas-Hesse 1998, private communication), 
which would correspond to a downward correction by a factor 
of — 1.7 for e.g. 14^(4686). In our observations the contribu- 
tion from an underlying population is, however, most certainly 
smaller since they are taken with smaller apertures including 
the bulk of the young emission line regions (cf. also Schaerer 
et al.l997) but only a smaller fraction of the continuum light. 

Given these potential difficulties both comparisons of rel- 
ative line intensities W-R/H/3 and W-R equivalent widths will 
be performed for all objects. Significant differences between 
the two methods likely indicate some difficulty with (d) and/or 
effects (i)-(iii). In this case the simplest meaningful compar- 
ison is between the observed WR equivalent widths and the 
maximum value predicted by the models (irrespectively of the 
age); it only depends on (Hi). 

It is well known that the W-R populations strongly depend 
on metallicity. We will therefore compare the observed regions 
to the models at the closest metallicity available. For the re- 
mainder of this Section we will associate the objects (O/H 
given in Sect, ^j]) with the following metallicity Z: NGC 3049: 
Z <; 0.02 (star symbol in Figs. [lO|to He 2-10: Z = 0.008 
(open triangle), and Z = 0.004 for the remaining objects (NGC 
3125, NGC 5253, Tol 89; all filled symbols). 

5.2. Results 

In Fig. |l^ we compare the observed relative W-R line intensi- 
ties of all W-R regions with the model predictions of SV98 at 
Z = 0.02 (solai-), 0.008, and 0.004 for an instantaneous burst 
with a Salpeter IMF (hereafter "standard" model; variations ai^e 



alent widths in Fig. |ll|. In both plots. Win/]) is used as a time 
indicator (decreasing W with time) in the synthesis models. 
The following general comments can be made to Figs. |l^ and 

0- 

From object to object the observed line ratios vary by a 
much larger factor than the W-R equivalent widths. This can 
be explained by the more rapid decrease of the Lyman contin- 
uum luminosity (i.e. the H/3 intensity) compared to the smaller 
variations in the optical continuum light in a short burst. With 
the exception of C IV A5808 in Tol 89 (cf. below), larger line 
ratios are found in objects with higher Z. This is to first or- 
der understood by the more rapid decline of the ionizing lumi- 
nosity with increasing metallicity, whereas the increase of the 
W-R population is less important (SV98). These findings fa- 
vor the picture of short burst periods (compared to the lifetime 
of massive stars) and support the predicted dependence of the 
ionizing pop ulation with metallicity fou nd independently from 
other studies ( Garcfa-Vargasetal., 1995 ; Stasinska & Leitherer, 
1996). 

W-R signatures are observed at H/3 equivalent widths 
where the models indeed predict W-R stars. Although this is not 
a strong constraint it indicates that the ly (H/3)-age chronome- 
ter is reasonably synchronised to first order Some difference 
seems, however, to be present for the appearance of WC signa- 
tures (see also Fig.jl^). Given the uncertainty (b) in the models, 
this difference is not significant (see below). 

Comparing the individual objects (classified by metallicity) 
with the model predictions, we find the following: 



NGC 3049 (Z ;> 0.02) 

The predictions of the standard model (instantaneous burst, 
Salpeter IMF) for the high metallicity region in NGC 3049 
show an excellent agreement with all observed W-R line 
strengths. This is also the case for the C III A5696 feature (not 
shown here) which is only detected in NGC 3049. 



He 2-10 (Z = 0.008) 

For He 2-10 the comparison between the equivalent widths and 
the line ratios shows a somewhat different picture. While the 
observed equivalent widths are considerably smaller than the 
maximum value predicted by the models (for the correspond- 
ing Z = 0.008), the relative line intensities are close to the 
predicted maximum. This difference can be attributed to the 
mismatch between gaseous and stellar emission (cf. (i)) seen 
in Fig. [l|, although other explanations (e.g. leakage of photons) 
cannot be ruled out. Assuming that effect ( Hi) is small (cf. be- 
low) we conclude from the comparison of the W-R equivalent 
widths (Fig. |ll]) that the W-R and O star populations are com- 
patible with an instantaneous burst and a Salpeter IMF. The 
same conclusion was akeady reached by Schaerer ( 1996 ) based 



Daniel Schaerer et al.: Populations of WC and WN stars in Wolf-Rayet galaxies 



15 




Fig. 10. Observed and predicted W-R line intensities relative 
to H/3 as a function of the H/3 equivalent width. The symbols 
used are show on the Figure. Filled symbols designate objects 
with a metallicity of Z ^ 0.004; open symbols have a larger 
Z (see text). The fairly small uncertainties on VF(H/3) are not 
shown in this Fig. (see Fig. |l^). Model predictions from SV98 
(instantaneous burst, Salpeter IMF) are shown for Z — 0.020 
(dotted), 0.008 (dashed) and 0.004 (solid). Upper right: relative 
line intensity for Hell A4686; lower left: CiV A5808; lower 
right: N III A4640+C III A4650 blend. Discussion in text 



NGC 3125, NGC 5253 and Tol 89 (Z = 0.004) 

Only some of the observed W-R line intensities and equiva- 
lent widths lie well on the predicted model curves plotted in 
Figs, p^and NGC 5253 represents the best case (Schaerer 
et al., 1997). For the considered objects we dispose in total of 
15 measurements of broad emission lines (all lines and all re- 
gions). The relative line intensity of 7 data points is above the 
maximum value predicted by the standard model. On the other 
hand the majority (12 of 15) of the equivalent widths are well in 
the range of the model predictions. The three remaining cases 
are W{He II A4686) of NGC 3125-A and Tol 89, which could 
well be overestimated due to nebular contamination and con- 
tributions from other lines in the blended region, and an excep- 
tionally large C IV A5808 in Tol 89 (cf. below). 

The uncertainties affecting the observed quantities have 
been discussed above (Sect. 5.1). Only for Tol 89 we have di- 
rect evidence that effect ( /) (displacement gas-stars) is likely of 
importance (Sect. ||); this does, however, not exclude that the 
line intensities in some other objects are also affected. Even af- 
ter reducing /(WR)/H/3 by a factor ^1.4 due to effect (//), 
discrepancies in the Une intensities remain. The observed W- 



Fig. 11. Same as Fig. |lO| for the W-R equivalent widths. Note 
that the observed value of W(C IV A5808) for Tol 89 exceeds 
the scale shown here (see Table ||). Discussion given in text 



though the observed values could be underestimated due to ef- 
fect We feel, however, that our measurements are only 
weakly affected by the old stellar populations (see Sect. 5.1). 



Other possibilities to explain the observations of some regions 
with large W-R/H/3 intensities and large VF(WR) may e.g. be 
a flatter IMF. These are discussed in Sect. |5.3| . 

Regarding the signatures essentially attributed to WC stars 
(i.e. N III A4640H- C III A4650, C IV A5808) we note a shift in 
W(Hf3) between the observed and predicted equivalent widths 
(bottom panels of Fig. |llj and Fig. [l^). The shift looks as if 
WC stars appeared to early and/or did not remain alive long 
enough. Possible explanations for this "chronometer-shift" are: 
effects (/), (b), (d), or variations in the burst duration, IMF, etc. 



(see Sect. 5.3). The uncertainties in the interpolation techniques 
determining the number of WC stars {(b)) seem large enough 
to explain this effect. 

For Tol 89 we measure an extraordinarily strong C IV 
A5808 emission (/(58G8)//(H/3) - 0.13, 1^(5808) ~ 12 A) 
which exceeds the values shown Figures [ll|and[l|. The other 
W-R emission lines do not show any particularity compared to 
the other regions at Z ^ 0.004. We have no simple explana- 
tion for this strong emission. Other observations are required 
to confirm the measurements for Tol 89. 



5.3. Age, duration and IMF of starbursts 

The mere detection of W-R signatures in the integrated spec- 
trum of galaxies reveals the presence of stars with ages gener- 
ally between 1 and 8 Myr (Meynet, 1995; SV98), depending 
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IMF; a = 2.35, 1. - Burst duration: 0, 1, 2, 4 Myr 




100 200 300 100 200 300 

W(HjS) W(H/S) 



Fig. 12. Observed and predicted W-R equivalent widths for ob- 
jects with Z ~ 0.004. Observations shown using the same sym- 
bols as in Fig. [l^. Models with different IMF slopes (left pan- 
els: a = 2.35 =Salpeter; right: a = 1.0) and burst durations of 
/Si — Q ("instantaneous", solid lines), 1 (long dashed), 2 (short 
dashed), and 4 Myr (dotted) are shown. Discussion given in 
text 



lutionary models, apply for stars which have undergone "nor- 
mal" mass loss. The possibility of additional mass loss through 
Roche-lobe overflow in a close binary system is not considered 
here ( Cerviiio et al., 1996|; [Schaerer & Vacca, 1996|; Van bev- 



eren et al., 1997; [Schaerer & Vacca, 1998| ; [Cerviiio, 1998| ). An 
independent age estimate is obtained from iy(H/3). From our 
standard models at the appropriate metallicity, we obtain ages 
between 3 and 6 Myr, which are reported in Table ^ (column 2) 
for each starburst region. As mentioned earlier these ages agree 
well with age range predicted for the W-R phase by the SV98 
models at the appropriate metallicity and give some support to 
V7(/3)-age relation. 

Models with different burst durations are plotted in Fig. [I^. 
These show that, within reasonable assumptions for the IMF 
slope, the duration of star formation episode is limited to A< ^ 
2 — 4 Myr. From Fig. [l^ and taking into account effect ( Hi) this 
is a conservative limit; longer durations are clearly excluded 
from this basis (see also Sect. 5.4). This quantitative estimate 



is compatible with most studies of W-R galaxies which have 
so far clearly fa voured instantaneous burst scenarios (Kunth 
& Sargent, 1981; [Arnault et al., 1989[ ; [Vacca & Conti, 1992[ ; 
Meynet, 1995] [Schaerer, 1996[ ). 

What can be said about the IMF slope? As mentioned ear- 
lier the observations of the W-R features in NGC 3049, He 2- 
10, and NGC 5253 ai^e quite compatible with a Salpeter IMF. 



A and Tol 89, and 14^(5808) in Tol 89 (representing in principle 
already lower limits due to ( Hi)) may require more W-R stars. 
This can be obtained by invoking a flatter IMF, which increases 
the WR equivalent widths as shown in Fig. |lj . However, we 
note that intrinsicly the model predictions for He II A4686 are 
not very sensitive to changes of the IMF slope. More IMF- 
sensitive are the features of WC stars, which descend on the 
average from more massive predecessors than WN stars and 
from a narrower range of initial masses. 

Although the standard models show some deficiencies it 
is difficult to claim significant differences with respect to a 
Salpeter IMF for the following reasons: 1) The Hell A4686 
predictions are not very sensitive to the IMF slope. Further- 
more the prediction of this line is uncertain (see (c), and the 
observed line may also be c onta minated by nebular emission 
(e.g. in NGC 3125-A, Sect. 2) Most of the WC signa- 
tures from the discrepant objects are within the range of the 
model predictions for a Salpeter IMF. We therefore conclude 
that within the uncertainties our observations are compatible 
with a Salpeter IMF. Although we cannot exclude this pos- 
sibility, no clear case requiring a significantly flatter IMF is 
found. Much steeper IMF slopes may, however, not be compat- 
ible with our data. 

Regarding the set of stellar tracks (cf. (a)) we note that 
the present results are only obtained for the stellar evolution 



models of Meynet et al. (1994) adopting high mass loss rates. 



Significantly lower equivalent widths of He II A4686 (approx- 
imately factor < 0.5), N III A4640H- Cm A4650 (< 1/6), and 
Civ A5808 (^ 0.1) are obtained with the standard mass loss 
models due to the reduced WN and WC populations (Meynet, 
1995). As mentioned before the high mass loss models are 
favoured by several independent comparisons with individual 
W-R stars and populations in the Local Group (Maeder & 
Meynet, 1994). This justifies their use when analysing massive 
star populations in extra-galactic objects. 



5.4. On the influence ofR136 Uke stars on spectra of W-R 
galaxies 

We shall now briefly digress and investigate how robust the 
above results (especially the burst duration) are in view of re- 
cent studies of massive stars in R136 and the Galactic Hll re- 
gion NGC 3603. observations of these regions show the 
presence of very young (^ 2 Myr) stars with considerably 
strong Hell A4686 emission (Drissen et al., 1995; de Koter et 
al., 1997). Although their exact spectral types are still under 
debate (e.g. Crowther & Dessart, 1998) quantitative analysis 
indicate that these stars are likely massive hydrogen-burning 
stars which are less evolved than typical WN stars (de Koter et 
al., 1997; Crowther & Dessart, 1998). According to the criteria 
used to define W-R stars, traditional synthesis models would 
most likely miss to count such stars as W-R and hence their 
He II A4686 emission would not be included. If such stars are 
common in young star forming regions, the He II emission pre- 
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Fig. 13. Predicted W-R equivalent widths (He II A4686) as a 
function of W^(H/3) for the standard instantaneous burst model 
(solid, SV98), and models including He II A4686 emission 
from massive (A^ini > 80 M©) main-sequence stars ("R136- 
like" stars) with varying burst durations: long-dashed (instan- 
taneous burst), dashed (At = 2 Myr), dotted (At = 4 Myr), 
dashed-dotted (continuous burst). All models are for Z=0.004 
and a Salpeter IMF. 



presence of such emission line stars reUeve the need to invoke 
very short star formation timescales ? 

To answer this question we investigated the effect of 
"R136-like" stars on the interpretation of W-R galaxy spec- 
tra, by performing test calculations (see Fig. To account to 
first order for He II A4686 emission from fairly unevolved stars 
(Olf and "R136-like" stars), we have assumed that all main- 
sequence stars with initial masses A/i„i > 80 M0 (cf. de Koter 
et al., 1997; Crowther & Dessart, 1998) show a He II emission 
with L4686 ~ 1-7 X 10'^^. This should provide an upper limit 
for the emission from "R136-like" starsQ. 

Figure [1^ shows that, fairly independently of the star for- 
mation history, the expected increase of He II A4686 due to 
the additional He II emitters is found only at T/F(H/3) ^ 200 
A. At lower equivalent widths the differences with the "stan- 
dard model" are quite small (cf. Fig. If a general phe- 
nomenon, the existence of "R136-like stars" may thus help to 
explain the observations of He II A4686 only in regions with 



^ SV98 included a contribution of Lmm = 2.5 x 10^^ ergs s~-^ 
from main-sequence stars with low gravity. For the three most lumi- 
nous stars of R136 and NGC 3603, the average 4686 luminosity is, 
however, larger: L4686 ~ (1-8 — 2.2) x 10^* ergs s^^ (Crowther, 
1997 and Drissen, 1997, private communications; de Koter 1996, pri- 
vate communication gives L4686 lower by approximately a factor of 



large H/3 equivalent widths. Our observations of W-R galaxies 
(see Fig. |lj and Table ^ and those of VC92 are mostly found 
at WiR/S) < 100 A. For these objects, the He II A4686 requires 
fairly short burst (see Fig. |lj and Fig. 2 of Schaerer 1996). In 
any case "R136-like" stars cannot explain the observed Nlll 
A4640H-C III A4650 and C IV A5808 emission which also re- 
quires short burst durations as discussed above. Other indepen- 
dent arguments against extended star formation in the observed 
W-R galaxies are discussed below. 

6. Summary and conclusions 

The main result of the present paper is the unambiguous de- 
tection of WC stars (indicated by broad C IV A5808 emission) 
in five previously known W-R galaxies, defined by broad He II 
A4686 emission mostly due to WN stars (Sect. 4.4). We con- 
firm the presence of WC stars in He 2-10 indicated by VC92. 
With our four new detections (in NGC 3125, NGC 5253, Tol 
89, and NGC 3049) the total number of extragalactic objects 
known to harbour both WN and WC stars is now 19 (cf. 
Schaerer & Contini 1998), which represent only ^ 20% of the 
total sample of W-R galaxies. The relative weakness of C IV 
A5808 compared to Hell A4686 (/(5808)//(4686) ^ 0.5 typ- 
ically) and its larger width requires sufficiently high S/N (^ 40) 
to be detected. As akeady pointed out by Schaerer et al.( 1997| ) 
this explains most likely the non-detection in previous obser- 
vations. 

In all objects broad lines of N III A4640H- C III A4650, He II 
A4686, and C IV A5808 are measured (Sect. A marginal 
detection of Hell A5412 is found in two W-R regions; weak 
C III A5696 indicative of late WC stars is found in the high 
metallicity W-R region of NGC 3049. From these emission 
lines we conclude that all W-R regions (except NGC 3049) con- 
tain a mixed population of WNL, and early WC and/or W03-4 
stars. This agrees well with expectations (e.g. Maeder, 1991; 
SV98). 

We have performed a detailed spatial analysis of the nebu- 
lar and W-R emission lines and the continuum light along the 
slit position in all our objects (Sect. H). In He 2-10 and Tol 89 
we found multiple peaks of nebular emission and a spatial off- 
set between the main peak and the stellar continuum. These 
structures are likely due to the existence of bubbles and loops 
in the ISM powered by the kinetic energy released by mas- 
sive stars (stellar winds and/or SNe, see Sect. 4.6). The spatial 
distribution of W-R stars follows closely the continuum and no 
significant distinction is found between WN and WC stars. The 
only exception is a bright He II A4686 peak with no continuum 
and nebular counterpart found in NGC 5253 (Schaerer et al., 
1997). Its origin is still unclear 

From the luminosity of the W-R signatures we have esti- 
mated the absolute number of W-R stars of the different sub- 



types (Sect. 4.5). For the regions whose H/3 luminosities vary 
approximately up to a factor of 20 the total W-R star content 
varies from ^ 30 to 1500. The estimated WCAVN number ra- 
tios (lower limits) are between 0.15 - 0.65, with typical val- 
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our objects with metallicities Z ^ 1/h Z©, these values are 
larger than the obse rved WCAVN ratios in L ocal Group ob- 
jects with similar Z ( Massey & Johnson, 1998 ). We argue that 
our WCAVN values are compatible with expectations for re- 
gions of short star formation. For He 2-10 and NGC 3049 the 
derived WCAV N ratio is below the trend given by Massey & 
Johnson ( 1998 ). This can have several explanations: 1) a short 
burst is observed at a particular time (quite low probability), 2) 
the number of WC stars is systematically underestimated (see 
Sect. \.5). The solution awaits new observations and quantita- 
tive analysis of "WCh-WN galaxies" of different metallicities. 

The detection of emission in the "traditional" blue W-R 
bump yields useful information about the age of the stellar pop- 
ulation and the presence of massive stars (cf. Sect. Above 
this, detecting both WN and WC features provides a consid- 
erable improvement for the following reason: WC stars are 
strongly evolved descendents of massive stars revealing He- 
burning products on their surface; prior to this phase these ob- 
jects are WN stars showing H-burning products (e.g. Maeder 
& Conti 1994). Given this WN WC sequence, which is ex- 
pected to be followed only by the most massive WN stars, it is 
clear that the predictions of WCAVN populations are particu- 
larly sensitive to the evolutionary scenario and burst parameters 
(e.g. IMF, burst duration). In Sect. || we have exploited this fact 
by performing detailed comparisons of three different observa- 
tional W-R signatures from WC and WN stars with predictions 
from the recent synthesis models of SV98. 

The comparisons of the observed equivalent widths and Une 
intensities relative to H/3 with models do not all show a sim- 
ple picture. The most important effects and uncertainties which 
may affect such a comparison have been amply discussed in 
Sect. H Some W-R signatures in few regions exceed the pre- 
dicted relative intensities and/or equivalent widths. The major- 
ity of the observed quantities (N III A4640H-C III A4650, He II 
A4686, C IV A5808, and H/3) can, however, be reproduced rea- 
sonably well by the SV98 models with a Salpeter IMF. Al- 
though some W-R lines may indicate a flatter IMF in some 
regions, no clear case requiring a significantly flatter IMF is 
found. Much steeper IMF slopes may, however, not be com- 
patible with our data. These results are in agreement with other 
studies of similar objects (e.g. Mas-Hesse & Kunth 1991, 1998; 
Schaerer 1996; Leitherer 1998). 

In order to reproduce the W-R lines, young populations 
with short durations of star formation are required. From our 
quantitative modeling, we find a conservative limit for burst 
durations of typically At 5s 2 — 4 Myr. A simple experiment 
shows that this result holds even if very young massive emis- 
sion line stars such as found in R 136 and NGC 3603 are com- 
mon in our observed star-forming regions. Such short star for- 
mation timescales can be understood if the light observed in 
our W-R regions comes from one or few individual compact 
regions, such as the super-star clusters frequently identified on 
high resolution //5r images (e.g. Conti & Vacca, 1994; Meurer 
et al., 1995). The finding of young and short bursts is supported 
by other independent constraints from Hll galaxies which have 



served spread of VF(H/3) and variations of nebular line ratios 
dStasiriska & Leitherer, 1996| ). 

Last, but not least, we also note from our comparison with 
starburst models that it is not possible to reproduce the ob- 
served WN and WC signatures adopting evolutionary models 
using standar d mas s loss rates. The high mass loss models of 
Meynet et al.( 1994 ) are clearly favoured from our comparison 
(see also Maeder & Meynet, 1994). A similar conclusion was 
obtained from comparisons of new evolutionary tracks with the 
observed W-R population in the extremely metal-poor galaxy 
I Zw 18 ( pe Mello et al., 1998| ), showing the usefulness of 
such observations to constrain evolutionary scenarios in envi- 
ronments (i.e. extreme metallicities) inaccessible in the Local 
Group. 

Our successful finding of WC stars in W-R galaxies opens 
the door to new systematic studies of massive star populations 
in starbursts. In addition to the subjects addressed in the present 
paper the study of WC stars in starburst regions is of interest for 
several other reasons: 1) WC stars may contribute to a harder 
ionizing spectrum which has implications on the observed neb- 
ular properties (e.g. Terlevich & Melnick 1985, Schaerer 1996). 
2 ) Regions which harbour WC stars during their lifetimes are 
expected to produce ejecta of significantly different compo- 
sition (cf. Maeder 1992). Our understanding of stellar evolu- 
tion, chemical evolution and the starburst phenomenon should 
greatly benefit from such future studies. 
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